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ABSTRACT 



Oh 



We investigate the properties of a clump-cluster galaxy at redshift 1.57. The morphology of this galaxy is dominated 
by eight star-forming clumps in optical observations, and has photometric properties typical of most clump-cluster and 
chain galaxies. Its complex asymmetrical morphology has led to the suggestion that this system is a group merger of 
several initially separate proto-galaxies. We performed Ha integral field spectroscopy of this system using SINFONI 
on VLT UT4. These observations reveal a large-scale velocity gradient throughout the system, but with large local 
kinematic disturbances. Using a numerical model of gas-rich disk fragmentation, we find that clump interactions and 
migration can account for the observed disturbed rotation. On the other hand, the global rotation would not be expected 
for a multiply merging system. We further find that this system follows the stellar mass vs. metallicity, star formation 
rate and size relations expected for a disk at this redshift, and exhibits a disk-like radial metallicity gradient, so that 
the scenario of internal disk fragmentation is the most likely one. A red and metallic central concentration appears to 
be a bulge in this proto-spiral clumpy galaxy. A chain galaxy at redshift 2.07 in the same field also shows disk-like 
rotation. Such systems are likely progenitors of the present-day bright spiral galaxies, forming their exponential disks 
through clump migration and disruption and fueling their bulges. Our present results show that disturbed morphologies 
and kinematics are not necessarily signs of galaxy mergers and interactions, and can instead result from the internal 
, evolution of primordial disks. 

Key words. Galaxies: formation - Galaxies: kinematics and dynamics - Galaxies: high-redshift - Galaxies: evolution - 
Galaxies: interaction 
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1. Introduction surface brightness galaxy (LSB) progenitors at high red- 
shift, but their clumps are one order of magnitude brighter 

Giant molecular clouds and star formation complexes in than the sta r-forming regions of LSBs (Smith et al. 2001) 

present-day spiral galaxies have masses that rarely exceed an( j the c l um p S are not only bright star-forming complexes 

10 of the disk mass, and the total mass in these low- but are m assive structures also visible in the near infrared 

mass star forming regions does not dominate the disk mass. (Elmegreen & Elmegreen 2005, hereafter EE05, see also 

At high redshift (z ~ 1 and above) the dumpiness and NICMOS observations by Dickinson 2000). "Clump clus- 

asymmetry of galaxies increases (Conselice 2003; Conselice ter » ga l a xies have the appearance of highly clumpy disks, 

et al. 2004; Daddi et al. 2004), which affects both the light Some are bulgeless, others have small red central concen- 

fraction that is in the clumps and the individual masses trations that resemble primordial bulges. Because the pho- 

of these clumps. Clumps at high redshift can be kpc-wide tometrical properties of their clumps are similar to that 

and as massive as 10 9 M (Elmegreen 2007, and references of chain galaxies, Elmegreen et al. (2004) suggested that 

therein). In particular, the so-called "chain galaxies" have chains are the e dge-on counterpart of such clumpy sys- 

striking morphologies dominated by alignments of massive temSj which solves the pro blem of the face-on counterpart 

clumps (Cowie et al. 1996; van den Bergh et al. 1996; of chains being m i ss i ng (O'Neil et al. 2000). The fraction 

Moustakas et al. 2004). Dalcanton & Shectman (1996) first of cnains and c i ump clusters increases with redshift, and 

proposed that these chain galaxies could be edge-on low because these galaxies do not have an exponential profile 

already established, Elmegreen et al. (2005) proposed that 

Send offprint requests to: F. Bournaud e-mail: they could be the progenitor of Z ~ 1 spirals that further 

frederic.bournaud@cea.fr evolved into present-day bright spirals. They would thus be 

* Based on observations obtained at the Very Large Telescope a crucial step in the formation of massive disk galaxies. 
(VLT) of the European Southern Observatory, Paranal, Chile 
(ESO program ID 278.A-5009) 
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The origin of the massive clumps in chains and clump 
cluster galaxies, however, remains uncertain and largely de- 
bated. Noguchi (1999) and Immeli et al. (2004a, b) proposed 
that these could be dense fragments resulting from Jeans 
instabilities in gas-rich and bulgclcss disks. Bournaud, 
Elmegreen, & Elmegreen (2007, hereafter BEE07) further 
explored this hypothesis and confirmed that these would 
evolve into spiral disks with a central bulge and an expo- 
nential disk (generally truncated in the form of type II as 
defined in Pohlen et al. (2007)). This hypothesis of internal 
instability however lacks clear support from observations 
so far. Spectroscopic observations of high-redshift disks 
unveil high turbulent speeds (Genzel et al. 2006; Forster 
Schreiber et al. 2006) so that the Jeans mass should be 
high and massive clumps could be formed, but these ex- 
tensive spectroscopic surveys lack associated deep imag- 
ing to directly adress the origin of clumpy galaxies. In 
contrast, Taniguchi & Shioya (2001) favored a multiple- 
merger origin involving several initially independent proto- 
galaxies. In that vein, clump-clusters could be the high red- 
shift equivalent of local compact groups, and some indeed 
have Stefan's Quintet-like morphology (at the resolution of 
high-reshift observations). The merging subunits could be 
aligned along filamentary structures in the case of chains 
(Taniguchi & Shioya 2001). Even binary major mergers can 
look clumpy with several star forming knots (Goldader et 
al. 2002), and even an exponential profile would not ex- 
clude a merger origin (Moran et al. 1999). Complex mor- 
phologies with severe disturbances and asymmetries, like 
for instance bent chains or shrimp-like systems (Elmegreen 
& Elmegreen 2006; Elmegreen et al. 2007a) further suggest 
that and interaction/merger origin is plausible at least for 
these categories. 

In this paper, we present kinematical IFU observations 
of the prototype clump-cluster galaxy UDF 6462 in the 
Hubble Ultra Deep Field (UDF) . The morphology and pho- 
tometry of this system at z = 1.57 have been studied in 
EE05 and Elmegreen et al. (2007b). A color image and 
an HST/ACS V-band image from the Hubble Ultra Deep 
Field (UDF, Beckwith et al. 2003) are presented in Fig. 1 
and 2, the main photometrical properties are recalled in 
Sect. 2.1. This system is one that a priori poses a problem 
for the hypothesis of clumps being internally formed by 
instabilities in a massive disk: it is globally asymmetrical 
and has a Southern extension that resembles a spiral arm. 
Furthermore, two nearby systems (in projection) could be 
involved in an on-going interaction. UDF 6462 is classified 
as an on-going merger by Conselice et al. (2007). 

Integral field spectroscopy data obtained with SINFONI 
on UT4 at ESO/VLT reveal complex internal kinematics, 
which has local disturbances but is overall dominated by a 
velocity gradient resembling that expected for an internally 
rotating system. Using the numerical simulations presented 
in BEE07, we propose a model of an unstable lopsided disk 
which forms internal clumps and, under some projections, 
a morphology resembling that of UDF 6462. We show that 
the observed velocity gradient and kinematical disturbances 
correspond precisely to those expected in such a scenario. 
Thus, this object that could a priori be regarded as being an 
on-going merger prototype, actually has properties that are 
fully compatible with the internally unstable disk hypothe- 
sis. Furthermore, the kinematics and some other properties 
are more difficult to explain in a multiple-merger scenario. 
Kinematical results of another nearby clumpy galaxy with 



a bent-chain morphology (UDF 6911) at z ~ 2 are also pre- 
sented, and also show a large-scale velocity gradient favour- 
ing an internal disk fragmentation rather than a merger 
origin. 

SINFONI observations of UDF 6462 are presented in 
Section 2. In Section 3, we present the numerical model 
that helps to interpret the morphological and kinematical 
properties of this system. In Section 4, we discuss the origin 
and evolution of this clump cluster galaxy, and our con- 
clusions are summarized in Section 5. Spatially resolved 
SEDs in this system are presented in Appendix A ; nu- 
merical models of metals distribution in disks and merg- 
ers are discussed in Appendix B. In this paper we use a 
Chabrier (2003) IMF and standard cosmological parame- 
ters with H = 75 km s _1 Mpc -1 . Throughout the paper, 
we will call clumps the large and dense structures typical 
of high-redshift clump-cluster galaxies, that have masses of 
10 8 " 9 M and typical sizes of the order of a kpc, much 
more massive than kpc-sized star- forming " clumps" of low- 
redshift spiral galaxies. 



2. Observations 

2.1. The UDF 6462 clump cluster 

UDF 6462 (Fig. 1) has been imaged in B, V, i and z 
bands using HST/ACS and J and H bands with NICMOS. 
The main properties are summarised in Table 1. A spec- 
troscopic redshift of z = 1.570 has been measured during 
the FORS2 campaign on the GOOD-S field (Vanzella et al. 
2006, 2008) based on detection of the [Oil] line. BViz data 
from Giavalisco et al. (2004) are presented and analyzed by 
EE05 and a photometrical analysis including NIR bands is 
in Appendix A. We here show V- and H-band images in 
Fig. 2 to recall the morphology of the system. EE05 esti- 
mate a stellar mass of 3.3 x 10 10 Mq and have identified 8 
star-forming clumps in this system, with 36% of the total 
2775 luminosity in these clumps and 64% in a diffuse and 
redder inter-clump medium. The typical age of the clumps 
derived from optical photometry is around 310 Myr, while 
the inter-clump medium has a typical stellar age of nearly 
3 Gyr. These properties, as well as the color-magnitude 
properties of the clumps, are globally typical of the sample 
of ten clump-cluster galaxies studied by EE05. The clump- 
cluster also has a central reddish blob, and the NIR emis- 
sion is much more concentrated than the optical emission; 
this resembles a primordial bulge. In the following we will 
call this structure "bulge" and subsequently in the paper we 
will make the case for this; for the moment we only imply 
this word to be indicative of a location within the system 
and not presuming of its real nature. 

The B — i ACS image in Fig. 3 shows a color gradient from 
this red central proto-bulgc to the bluer clumps. This is not 
simply reddening caused by higher extinction, considering 
that the region with red B — i colors corresponds to the 
area of highest H-band surface brightness (see also SED 
fits in Appendix A). Photometric measurements compared 
with stellar population models suggest a bulge mass frac- 
tion of about 15-20%. This is actually an upper limit be- 
cause of the NICMOS H-band resolution being larger than 
the red blob diameter. Spectral energy distributions of this 
central bulge-like blob and the other clumps are shown in 
Appendix A and also discussed in Section 4. 
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Fig. 2. The UDF 6462 clump-cluster: (a) HST-ACS V-band image plotted in log-scale. UDF 6911 and 6886 are are very different 
redshifts (see text) and not interacting with UDF 6462. Dashed arrows indicate the pseudo-slits used to derive position-velocity 
diagrams, (b) NICMOS H-band image plotted in log-scale, showing a more concentrated emission around the central bulge-like 
blob, (c) H-band continuum from our SINFONI observations plotted in log-scale: UDF 6462 is detected, including its faint southern 
extensions; UDF 6911 is also detected. 
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Fig. 1. Color image of the UDF 6462 clump-cluster from optical 
multi-band ACS imaging. The bent-chain UDF 6911 is visible 
to the left. 



UDF 6462 has two neighboring systems, which we here 
denote UDF 6911 and UDF 6886 (see labelling in Fig. 2). 
UDF 6911 has a chain- like morphology, while UDF 6886 
is more compact. The colors and dimensions are compa- 
rable to that of UDF 6462, possibly suggesting a similar 
redshift. We actually find from our SINFONI observations 
(see below) that these are at different redshifts. 

The star formation rate (SFR) derived from the 
Spitzer MIPS 24 /um flux (Chary et al. in preparation) is 
85 Mq yr _1 based on the Chary & Elbaz (2001) models. 
This value could be an over-estimate because the 24 /xm 
flux is likely contaminated by nearby systems. The UV- 
based estimate corrected for extinction is 50 Mq yr -1 , 
as in Daddi et al. (2007), and is not contaminated by 
nearby systems. The [Olljflux, corrected for extinction, cor- 
responds to an SFR of 35-55 Mq yr _1 (using relations 



UDF 6462 



a (J2000) 


3h 32m 43.65s 


5 (J2000) 


-27° 46' 59.2" 


redshift 


1.571 


B mag 


24.16 


i mag 


23.60 


M» 


3.3 x 10 10 M 


SFR,uv 


50 M Q yr" 1 


SFRir 


85 M yr" 1 


12+log(0/H) 


8.53 ±0.02 


Vcirc 


~ 100 km s" 1 



Table 1. Main properties of UDF 6462. Details and references 
are in Section 2.1. 
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Fig. 3. B - i color map of UDF 6462 from HST/ACS data, 
showing the color gradient from the central proto-bulge to the 
outer clumpy regions. 



from Kennicutt 1998a). As the 24 /im flux may be con- 
taminated by nearby systems (in particular UDF 6911 and 
6886), these various estimates are consistent and suggest 
an SFR of about 50 Mq yr -1 with an uncertainty around 
20 Mq y^ 1 . 
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2.2. SINFONI observations and data reduction 

The observations were carried out using the SINFONI near- 
infrared integral field spectrometer mounted on the VLT 
UT4 telescope (Eisenhauer et al. 2003; Bonnet et al. 2004; 
Gillcssen et al. 2005). An integration time of 5 hours was 
spent on the target in December 2006 with average seeing 
of 0.5 to 0.6" in the H-band. No adaptative optics was used 
because no suitable Natural Guide Star was found close 
enough to the target. We used the high-resolution H-band 
grism, which offers a FWHM spectral resolution of R ~ 
3000 in the H-band, equivalent to a FWHM of ~ 90 km s" 1 
for the Ha line at z = 1.57. A first standard reduction of 
the data has been performed using the SINFONI pipeline 
(Modigliani et al. 2007). 

We then performed a more careful reduction relying on a 
set of standard IRAF tools to reduce longslit spectra (Tody 
1993), which we modified and extended by a set of IDL rou- 
tines matched to the special requirements of SINFONI. The 
dark-frame subtracted and flat-fielded data are rectified us- 
ing an arc lamp (for the spectral direction) and an artificial 
point source (for the spatial direction) . To account for some 
spectral flexure between frames, we rectify the data before 
sky subtraction. We correct for variations of the sky back- 
ground in the object and the sky frame by rescaling the 
count rate of the sky frame to that of the corresponding 
object frame, carefully masking the object, before recon- 
structing, aligning, and combining the three-dimensional 
data cubes. Flux scales are obtained from standard star 
observations taken once per hour at an air mass and posi- 
tion which match those of the source. We also used stars to 
measure the spatial resolution. The spectral resolution was 
measured from night sky lines. For a more detailed descrip- 
tion see, e.g., Nesvadba et al. (2006). The data presented 
hereafter are those obtained with this careful reduction, but 
we checked that all features discussed in the paper are also 
seen in the dataset reduced with the standard SINFONI 
pipeline procedure. 

2.3. Spectroscopic redshift and neighbouring systems 

We identified the Ha, [Nil] and [SII] emission lines associ- 
ated with UDF 6462 in the SINFONI datacube; the associ- 
ated spectroscopic redshift is z = 1.571. In UDF 6911, we 
detect the [OIII] and H/3 emission; the associated redshift 
is z = 2.071. As for UDF 6886, a single emission line is 
detected and assumed to be the Ha line; in this case the 
redshift is z = 1.427. Results on the internal kinematics of 
UDF 6911 are shown in Section 4.3. The emission line of 
UDF 6886 is not spatially resolved in our data. 

The integrated spectrum of UDF 6462 is shown in Fig. 4. 
The HII emission line from UDF 6462 falls far enough from 
telluric OH lines so that robust moment maps can be ex- 
tracted. The detection of [Nil] and [SII] emission is fainter 
and limited to the central regions. Because these lines are 
faint and/or contaminated by OH lines, associated velocity 
fields could not be properly retrieved over a significantly 
extended area. 

The H-band continuum map of the system is shown in 
the third panel of Fig. 2. This continuum map is based only 
of the parts of the spectrum that are not affected by an OH 
sky line or another emission line. Linear interpolations were 
used to recover the continuum over the regions affected by 
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Fig. 4. Integrated spectrum of UDF 6462 central regions (top). 
Sky lines are shown below. Shaded bands corresponds to wave- 
length ranges where the sky emission level exceeeds the av- 
erage 1-cr noise level estimated on the continuum. The HII 
line from UDF 6462 falls in a region free of telluric OH lines 
for at least ±600 km s" 1 around the central wavelength. The 
[NII]a6582 line is also unaffected by sky lines. The 2-D spectrum 
shown on the middle corresponds to a pseudo-slit ranging from 
S = -27° 47' 58.5" (top) to 8 = -27° 47' 59.5" (bottom); the 
pixel size on this plot is the raw SINFONI pixel size in both the 
spatial and spectral directions. 

these lines, and the continuum intensity was finally com- 
puted for wavelengths within the H-band range. 

The two companions of the main system have largely 
different redshifts. In spite of their apparent proximity, sim- 
ilar colors and comparable clump sizes, these are not nearby 
systems in the 3-D Universe. The spurious apparent prox- 
imity and similarity of these systems illustrates that the 
identification of interacting/merging galaxies at high red- 
shift without accurate spectroscopic redshifts cannot be rig- 
orous. In the present case, the disturbed morphology of 
UDF 6462 could have further suggested an on-going inter- 
action with its apparent companions, which spectroscopic 
redshifts rule out, showing that morphology is not always 
a good tracer of on-going galaxy interactions. 

UDF 6462 is not interacting with its apparent neigh- 
bors; nevertheless this does not rule out a merger origin 
for this system itself. The blobs gathered into this irregular 
clump cluster may still be several independent units merg- 
ing together instead of internal fragment of a single system, 
and this is what is studied in the following. 

2.4. Results: kinematics of the UDF 6462 clump cluster 

The first moment map associated with the Ha emission 
line in UDF 6462 is shown together with the Ha intensity 
contours in Fig. 5. The main regions are detected with a 
signal-to-noise 1 (S/N) ratio above 3. The detection of the 
Southern blob is more uncertain, at about S/N=2.6, but 
this blob is also seen on the Position- Velocity (P-V) dia- 

1 we measured the S/N ratio as the ratio of the flux in the 
HII line to the noise measured on the continuum, following for 
instance Flores et al. (2006) with a single line in our case. 
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gram shown in Fig. 6, on the continuum map (Fig. 2). The 
velocity field in the main region shows a large-scale veloc- 
ity gradient with the Northern parts approaching and the 
Southern parts receding. The (suspected) Southern blob 
follows the same trend. One can however note that this 
velocity gradient is severely disturbed compared to what 
would be expected for a usual rotating disk. The V/a ratio 2 
is typically between 1 and 2. In particular, Fig. 5 shows sev- 
eral clumps with velocities larger or smaller than expected 
for a purely rotating system. The P-V diagram along the 
apparent major axis of UDF 6462, shown in Fig. 6, confirms 
the overall trend for a velocity gradient throughout the sys- 
tem but with local disturbances; in particular a clump at 
(5=-27°47'59.4" in the observation on Fig 6 has a velocity 
lower than expected from the large-scale gradient, and a 
similar perturbation is found in the model at X=4 kpc in 
the same figure. 

The kinematics of UDF 6462 thus shows some sign of 
rotation, with on average a velocity gradient along the 
major axis of the system. The velocity dispersion map 
shown in Fig. 6 has a maximum in the dispersion around 
6 ~ — 27°46'58.8" which is the region where the velocity 
gradient on the P-V diagram is maximal. This is a sig- 
nature expected for rotation (e.g., Forster Schreiber et al. 
2006; Genzel et al. 2006; Puech et al. 2006), in particular 
because the high-dispersion region is extended in the direc- 
tion perpendicular to the velocity gradient (i.e., along the 
morphological minor axis). Moreover, the highest slope in 
the P-V diagram lies on the continuum peak, which further 
indicates that this velocity gradient traces a rotation curve 
rather than some other kind of motion. 

This underlying large-scale rotation pattern is, however, 
largely disturbed at the local scale of individual clumps. 
Some have velocities smaller or larger than expected for 
pure rotation, with discrepancies up to 40-50 km s _1 . This 
is larger than the velocity anomalies caused by spiral arms 
and other internal structures in classical low-redshift disk 
galaxies (e.g., Rots 1975). As a result, the velocity disper- 
sions and local kinematical disturbances are high, nearly 
comparable to the rotation velocity if the large-scale gradi- 
ent is tracing rotation. 

2.5. Metallicity and central enrichment in UDF 6462 

Intensity maps of the [Nil] and [SII] emission lines are dis- 
played in Fig. 7. The [Nil] and [SII] emissions seem more 
centrally concentrated than the Ha emission. Metallicity 
estimates, based on [Nil] /Ha, within an aperture centered 
on the surface brightness peak in the H-band compared to 
that of the entire system, suggest that this region of high 
metallicity is real and not due to lower signal to noise in 
the outer regions of emission. 

We used the N2 ratio, defined by N2 = 
log (NII,\e584/Ha), as a proxy for the Oxygen abun- 
dance, as proposed by Storchi-Bcrgmann et al. (1994). 
Pettini & Pagel (2004) determined the relation: 

12 + log(0/H) = 8.90 + 0.57 x N2 (1) 

The N2 ratio in UDF 6462 supports star formation as 
the powering source of the optical lines (Erb et al. 2006a; 

2 The V/a ratio here is not corrected for inclination because 
the inclination is unknown (or model-dependent), so it is not a 
direct indicator of the rotational support in this system. 



Halliday et al. 2008), consistent with the fact that this sys- 
tem does not show a mid- infrared excess typical of AGNs. 
The N2 ratio is also not as high as in could be for sys- 
tems with strong winds (Lehnert & Heckman 1996). To 
derive the bulge metallicity, we measured intensities inside 
a 0.6"diameter aperture centered on the NICMOS H-band 
emission peak. The global metallicity was derived inside the 
3-ct contour of the Ha flux map. We find: 

- 12 + log(0/H) = 8.59 ± 0.04 for the bulge 

- 12 + log(0/H) = 8.53 ± 0.02 for the whole system 

where the indicated uncertainties correspond to the esti- 
mated noise level in emission maps 3 . This corresponds to a 
metallicity of Z ~ 0.5Zq for the whole system. The cali- 
bration of the N2 ratio might suffer from systematic uncer- 
tainties, but there is rather robust evidence for central en- 
richment in this clump-cluster system. The relative Oxygen 
abundance in the central redder blob is about 20% higher 
than in the rest of the system. The whole system agrees 
with the mass-metallicity relation of z ~ 2 established by 
Erb et al. (2006a) and Halliday et al. (2008) and has a 
slightly lower metallicity than Luminous Infrared Galaxies 
at z ~ 0.6 (Liang et al. 2004), but the bulge metallicity 
alone is significantly enriched and would be in between the 
z ~ 2 relation and that of low redshift galaxies (SDSS, 
Trcmonti et al. 2004). The metallicity map obtained from 
the N2 ratio displayed in Fig. 7 shows this central enriche- 
mcnt, and a radial metallicity gradient spanning from the 
central red bulge towards the rest of the system. 



3. Numerical Modelling 

Large kinematical disturbances are frequently interpreted 
as evidences for interaction/mergers, in particular when 
they are associated with complex morphologies (see 
Introduction). In the following, we use an V-body numeri- 
cal model to study whether an internal origin for the star- 
forming clumps in UDF 6462 could be another viable ex- 
planation of such properties. 

3.1. Simulation technique 

The numerical technique used here is similar to that pre- 
sented in more detail in BEE07 and references therein: we 
here only recall the main aspects. The N-body code models 
stellar, gaseous, and dark matter particles. The gravita- 
tional potential is computed on a 256 3 grid using an algo- 
rithm based on Fast Fourier Transforms and optimized for 
vector computers. The grid cell size and gravitational soft- 
ening length are 180 pc. We use 2 x 10 5 stellar particles, 
2 x 10 5 gas particles, and 3 x 10 5 dark matter particles. Star 
formation is accounted for using a local Schmidt law with 
an exponent 1.5 (Kcnnicutt 1998b). This 1.5 exponent is 
applied to the gas volume density in our models, which is 
equivalent to applying it to the surface density assuming a 
constant thickness, as is the case in our initial conditions. 
The dissipative nature of the turbulent interstellar medium 
is modeled with a sticky-particle scheme. 



Because our data are not flux-calibrated, these error bars 
are only related to the spatial variations of O/H; a constant 
additional offset may affect both measurements in the same way. 
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Fig. 5. Moment maps for the Ha line from 
our VLT/SINFONI data, (a) m = in- 
tensity map, (b) m = 1 velocity field, 
and (c) m = 2 velocity dispersion field 
for UDF 6462. The 2-cr and 5-cr contours 
of the zero moment map are shown; the 
southern blob is detected at about 2.6-cr, 
with a continuum signal at the same lo- 
cation (Fig. 2). Instrumental broadening 
was subtracted from the dispersion map. 
UDF 6462 shows signs of rotation with 
a velocity gradient and a central peak in 
the velocity dispersion; however kinemati- 
cal disturbances are larger than in classical 
low-redshift disk galaxies. 
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Fig. 6. Left: Position- Velocity diagram of UDF 6462 along its apparent major axis. The thick pseudo-slit used to compute this 
diagram is shown in Fig. 2, data have been resampled to a twice smaller pixel size on this P-V diagram. Note the global velocity 
gradient corresponding to the velocity field in Fig. 5, the steepest parts of the gradient around the emission peak corresponding to 
the peak in the line-of-sight velocity dispersions. Right: Position- Velocity diagram for the model in Section 3, along its projected 
major axis, at full resolution and degraded to the spatial seeing, pixel size and instrumental velocity resolution. A global velocity 
gradient as in UDF 6462 is reproduced, together with large local disturbances and rather high velocity dispersions. Each clump 
in this modeled clump-cluster has an internal velocity gradient associated with its own spin, which could not be resolved at the 
SINFONI resolution. 
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Fig. 7. (a,b) [Nil] and [SII] intensity 
maps of UDF 6462. Metals appear to be 
centrally concentrated around the reddish 
proto-bulge and UDF 6462 shows sign 
of central enrichment illustrated by the 
metallicity map (c) where the O/H ratio 
is estimated from the [NIIJ/HII flux ratio 
(see text) . 



3.2. Model parameters 

The physical parameters of the model presented here are 
close to that of model 6 in BEE07, except that a mod- 
erate initial asymmetry is introduced in the system. The 
baryonic disk initially has a flat radial profile (uniform sur- 
face density) . This is because clump-cluster galaxies gener- 
ally do not have a spiral-like exponential profile, but gen- 
erally show irregular and less concentrated profiles (e.g., 
Elmegreen et al. 2005), which is the case for UDF 6462 at 
least in optical bands. The hypothesis of a uniform surface 
density was made because it is a simple and not restrictive 
choice already made in BEE07 models, but is not required 
to produce a clump-cluster morphology - clump interac- 
tions and migration will anyway rapidly redistribute the 
disk masss. 



This disk has a radius of 8 kpc radius and a mass of 
6 x 10 10 Mq , with a constant scale- height of 700 pc. The ini- 
tial gas fraction is 50%, giving a stellar mass of 3 x 10 10 Mq 
4 . Stars are initially given random motions with a Toomre 
parameter Q s = 1.3, and gas particles have a velocity dis- 
persion of 11 km s _1 . 

As the central blob of UDF 6462 is redder than the rest 
of the system and could be a primordial bulge, a spheri- 
cal bulge of 14% of the disk mass is added to the initial 
conditions in the model. This bulge fraction is compatible 
with observational estimates for UDF 6462 (15-20% from 
the SED fits in Appendix A), because 14% is the initial 
fraction in the model, which grows with time during the 
clump cluster evolution, and because the observed value is 

4 This is compatible with the observational estimate (see 
Sect. 2.1) once some more stars have formed in the simulation 
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an upper limit to the actual bulge mass. Note furthermore 
that the bulge mass is not a crucial parameter to study the 
large-scale kinematics of clump cluster formed from frag- 
mented disks: even bulgeless models in BEE07 show the 
same kincmatical properties (disturbed rotation) as we find 
below for the present model. 

The dark matter halo is modeled as a Plummer sphere 
of total mass 5 2.9 x 10 11 Mq and scale-length 12 kpc, trun- 
cated at 24 kpc. The dark-to-visible mass ratio within the 
initial disk radius is 0.35. Models of unstable disk evolution 
with NFW dark matter profiles are presented in Elmcgrccn 
et al. (2008b), and the different halo profile does not rad- 
ically change the evolution of clump clusters formed this 
way. 

Because the observed system is asymmetrical, the disk 
center is spatially offset from the halo and bulge mass cen- 
ters by 10% of this radius. Given that the degree of mor- 
phological lopsidedness in present-day spirals is frequently 
as large as 10 or even 20% (e.g., Bournaud et al. 2005; 
Angiras et al. 2007; Reichard et al. 2007) this is a realistic, 
moderate initial asymmetry 6 . The first snapshot in Fig. 8 
illustrates the initial asymmetry. 

3.3. Morphology 

The gas-rich, unstable disk in the present model is similar 
to the one studied in BEE07 and follows a similar evolution. 
The disk fragments into massive clumps in about 10 8 years, 
the formation of which can be followed in Fig. 8 and 9. 
This system forms eight big clumps embedded into a fainter 
underlying disk; the mass fraction in the clumps reaches a 
peak value after t = 270 Myr when 41% of the disk mass is 
gathered in the clumps. 

Because of the initial asymmetry, the clump distribu- 
tion itself can appear quite asymmetric. Depending on the 
projections, the dumpiness can increase the apparent asym- 
metry. For instance, the last instant in Fig. 8 would have a 
high apparent asymmetry seen from the bottom-right, and 
a lower asymmetry seen from the bottom-left. We also note 
the development of an lopsided spiral structure from the 
initial asymmetry and instability. This structure resembles 
the Southern extension of UDF 6462, which might thus be 
formed without an external tidal interaction. 

The instability of this massive gas-rich disk, combined 
with its moderate initial asymtery, results in the quite 
asymmetrical and clumpy projection shown in Fig. 9 with 
an inclination of 64 degrees. This projection is morphologi- 
cally similar to the HST V-band image of UDF 6462 shown 
in Fig. 2. 

3.4. Kinematics 

We show in Fig. 9 the velocity field of the modeled clump 
cluster for the projection which looks similar to UDF 6462. 
The velocity field shows a large-scale velocity gradient, but 
local variations and irregularities can be seen which re- 
semble those observed in the SINFONI velocity field of 

5 This total mass is that of the un-truncated sphere. 

6 There are furthermore reasons to assume that z ~ 1 disks 
are more lopsided than z ~ 0. Be these asymmetries caused 
by interactions, by asymmetrical gas infall, or by halo asymme- 
tries, these causes themselves should be more frequent/efficient 
at high redshift. 



UDF 6462. In addition, the P-V diagram along the appar- 
ent major axis is shown in Fig. 6, both at full resolution 
and at the SINFONI resolution. Clumps on this P-V dia- 
gram show an internal spin that could not be detected at 
the SINFONI resolution. On larger scales, clumps follow a 
global velocity gradient associated with the rotation of the 
system, but the individual velocity of each clump can be 
offset from this global gradient by a few tens of km s^ 1 . 
Similar kincmatical irregularities are also shown for some 
cases in BEE07, and are caused by clumps interacting to- 
gether and migrating in the system. In the present case, 
the P-V diagram of the system also incorporates the asym- 
metry, which gives it an overall aspect qualitatively similar 
to the observed one shown in Fig. 6. A velocity gradient 
can be detected and traces the overall rotation of the sys- 
tem, but the internal evolution of the clump cluster causes 
significant disturbances to the kinematics, that are larger 
than the disturbances caused by spiral arms and internal 
structure in regular spiral disks (e.g., Rots 1975). 

The observed peak-to-peak velocity extension of the P- 
V diagram of UDF 6462 is about 120 km s _1 . Assuming 
that the circular velocity is traced by the envelope at 50% 
of the peak level (e.g., Donley et al. 2006), the projected 
circular velocity V c sini ~ 90 km s _1 . Assuming an inclina- 
tion i = 63 deg, as suggested by our numerical model and 
consistent with the observed morphology, this corresponds 
to a circular velocity of about 100 km s _1 . 

Both the morphological and kinematical properties of 
UDF 6462 are qualitatively reproduced by the clump- 
cluster model presented here, at the instant and under 
the projection that we have chosen. We do not claim that 
the real system has followed exactly the same evolution 
in particular, because we simply picked a similar-looking 
model among several ones, and the solution is probably 
not unique. Nevertheless, the present model illustrates the 
fact that the clumpy morphology and the disturbed kine- 
matics of UDF 6462 can be reproduced by models based on 
the internal evolution of unstable gas-rich disks: this model 
shows that the morphology, clumpy and asymmetric, and 
the kinematics with local perturbations, can both be ac- 
counted for by the internal evolution of gas-rich disks. To 
achieve this, we have here added a moderate initial asym- 
metry in the present model, which under a chosen projec- 
tion reproduces the peculiar morphology observed for this 
system and its disturbed morphology at the same time. 
Other models presented in BEE07 show similar morpholo- 
gies and kinematical disturbances of the same order, even 
if not reproducing the individual case of UDF 6462. 



4. Discussion 

4.1. UDF 6462: Unstable disk or multiple merger? 

The UDF 6462 clump cluster has been considered a good 
candidate for a merging system, because of its irregular and 
asymmetrical morphology compared to some other chain 
galaxies that have more symmetrical aspects. In such a 
group model, initially separate dark matter haloes with 
their stellar and gaseous content each would be merging 
together, and be presently observed in a phase resembling 
a compact group. It could also be a binary merger with 
several star forming clumps. 
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Fig. 8. Snapshots of the disk mass density in the model, with time indicated in Myr. The halo mass center is marked by the black 
+ symbol. The initial conditions are closely similar to the first snapshot, with a uniform disk density. The last snapshot shown 
here is slightly before the projection reproducing the UDF 6462 morphology (Fig. 9). Note that the asymmetrical aspect resulting 
from the initial disk-halo offset is amplified by a long spiral arm extension and the outer location of one of the clumps. 
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Fig. 9. Projection of the model repro- 
ducing the UDF 6462 clump-cluster mor- 
phology, after 380 Myrs of evolution, (a) 
projected density map of stars formed 
less than 300 Myrs ago to mimic the V- 
band map shown in Fig. 2. (b) line-of- 
sight velocity map. The velocity field was 
smoothed to a 1 kpc resolution (Gaussian 
smoothing) to reduce the numerical noise, 
and a cut-off along an isodensity con- 
tour was applied, (c) line-of-sight veloc- 
ity dispersion map. Under this projec- 
tion, model qualitatively reproduces the 
severely disturbed rotation pattern ob- 
served in UDF 6462, with a large-scale ve- 
locity gradient and a central peak in the 
velocity dispersion. (d,e,f) Same moment 
maps degraded to the spatial resolution of 
our SINFONI observations. Irregularities 
in the velocity field are still visible, for 
instance the high-velocity region around 
X=-2 and Y=-l kpc, which is compara- 
ble to the velocity anomalies observed in 
UDF 6462. 



However, our numerical model shows that such a com- 
plex morphology can result from the internal evolution 
of an unstable gas-rich disk galaxy. In this fragmenta- 
tion model, a single dark matter halo contains or accretes 
large amounts of gas, so that the primordial disk can be- 
come Jeans-unstable and fragments into several large star- 
forming clumps, taking the appearance of a clump-cluster. 
Only a moderate (realistic) initial asymmetry is needed in 
the model to qualitatively reproduce the projected mor- 
phology of UDF 6462. The kinematic data obtained with 
SINFONI unveil a large-scale velocity gradient, together 
with local disturbances. Such a large-scale velocity gradi- 
ent is a signature expected for a clump cluster formed by 
internal instabilities in a disk galaxy. The irregularities in 
the velocity field are then explained by clump-clump inter- 
actions that cause the individual velocity of each clump to 
differ significantly from the initial rotation velocity. Both 
the morphology and the kinematics of UDF 6462 can thus 
be accounted for by an unstable disk scenario. 



The viability of this first scenario does not strictly rule 
out a merger origin. The blobs observed in UFD 6462 could 
be initially independent sub-units merging together at the 
same time to form a more massive galaxy. However, several 
observed properties appear unlikely under this hypothesis: 



— (1) First, the clumps have comparable masses and 
sizes (see Figs. 2 and A.l) which is typical of most 
chains and clump cluster galaxies. The dynamical fric- 
tion timescale for merging galaxies is inversely propor- 
tional to the mass. Blobs that have the same mass as the 
most massive one should merge rapidly, and low-mass 
blobs should survive longer. The present morphology of 
UDF 6462 could then only be a transient one, unlikely 
to be observed. A massive blob surrounded by several 
lower-mass ones would be a more likely configuration. 
On the other hand, the similar masses for all clumps are 
naturally expected for the internal fragmentation of a 
primordial disk, because the typical clump mass in this 
case would be the Jeans mass in the outer disk. 
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— (2) The large-scale rotation-like velocity gradient and 
the central peak in the velocity dispersion are also char- 
acteristics of rotating disks that are predicted to be pre- 
served (with disturbances) in a clump cluster formed 
from an unstable disk. Such properties are not necessar- 
ily expected for a merging system. In particular, there 
is no reason for which merging blobs would be glob- 
ally rotating in the same direction. The observed ve- 
locity dispersions are rather high but do not dominate 
the large-scale velocity gradient; for a multiple merging 
they would more likely be larger the rotation velocity. 
Actually the merging timescale is faster for systems that 
rotate in the same direction, because of angular mo- 
mentum removal in tidal tails, while counter-rotating 
systems take longer to merge. The absence of counter- 
rotating unit among the several blobs is a possible but 
unlikely property for a multiple merger origin. 

— (3) The radial metallicity gradient (Fig. 7) with cen- 
tral enrichment in the reddish bulge-like blob is also 
expected for a (clumpy) disk, and the total metallic- 
ity is about that expected for a disk of this mass. The 
centrally concentrated [Nil] /Ha ratio is similar to what 
Forster Schreiber et al. (2006) found in a rotating disk 
at z ~ 2 (Q2343-BX610), with a ratio between the 
central metallicity peak and the outer disk metallicity 
larger by 35% but comparable to what we find here 
in UDF 6462. Such a disk-like metallicity distribution 
would be another coincidence if UDF 6462 is an-going 
group merger. Generic models of galaxy mergers and 
disk fragmentation in Appendix B confirm that a radial 
metallicity gradient is better preserved by the disk frag- 
mentation mechanism than in early-stage galaxy merg- 
ers with clumpy morphologies. 

— (4) Stellar population studies confirm that the clumps 
are young, star-forming objects around an older, central, 
bulge-like object. The observed SEDs and theoretical 
fits for the whole galaxy, the bulge, and several clumps 
are shown in the Appendix A and Fig. A.l. The central 
region that we have called bulge is clearly redder than 
the clumps. The reddest clump is an extended region 
jutting off from the bulge and may be part of the bulge. 
The central reddish bulge has an age around 1 Gyr. The 
clumps have an age of ~ 50 Myr with a comparable de- 
cay time in the star formation rate. There appears to be 
no significant old stellar component in the clumps other 
than the central bulge: if these clumps were separate 
galaxies in a merging process, they should have older 
stellar components in them, formed before this merger 
as in any galaxy. The absence of old stars in the clumps 
further supports our model in which they formed re- 
cently by the fragmentation of a gas-rich disk. Moreover, 
these clumps appear to have similar ages (similar SED 
slopes) which again arises naturally in the disk fragmen- 
tation scenario. 

The observed properties of UDF 6462 are thus consis- 
tent with an unstable-disk origin, while they would be un- 
likely coincidences in a merger origin. We also note that 
a stellar mass of 3.3 x 10 10 Mq and an external diame- 
ter of ~ 10 kpc (without the low-mass Southern extension) 
would place the system in agreement with the mass-size re- 
lation of disk galaxies from Trujillo et al. (2006) which, at 
a slightly higher redshift of 1.7, indicates an half-light di- 
ameter of 6 kpc for a galaxy of this mass. Furthermore, the 



5 kpc radius and 100 km s _1 circular velocity inferred above 
are consistent with the size-velocity relation observed for 
star-forming disks at z— 1.5-2 by Bouche et al. (2007), while 
submillimeter galaxies (likely major mergers) are found to 
have much lower size: velocity ratios by these authors. 

These various properties cannot definitely rule out a 
merger origin, also because only one system with direct ev- 
idence for massive clumps from HST/ACS and NICMOS 
imaging has been studied so far. However, for the merging 
of a group of initially separate galaxies, one would expect 
a concentration of dark matter around each visible blob, 
so that the largest velocity gradients should be internal to 
the clumps, and not in the form of a global rotation of 
the system. One might assume that, in a merger model, 
the initially separate dark haloes could have merged into 
a single one, while their baryonic counterparts would still 
be separate. However such a system should resemble Local 
Compact Groups, in which the prominent kinematical sig- 
natures are internal gradients inside each galaxy with lit- 
tle global rotation of the group (e.g., Amram et al. 2007, 
and references therein). Local systems have lower gas frac- 
tions, but tidal interactions in the early stages of mergers 
are dominated by the DM and stellar potential, so that the 
kinematic response of the gas would at first order be similar 
in higher redshift compact groups. This comparison further 
confirms that the properties of UDF 6462 arc unlikely for 
a merger mechanism. On the other hand, the fragmenta- 
tion model within a single dark matter halo predicts (see 
the simulation here and those in BEE07) that clumps have 
some internal rotation, but with low velocities (a few tens 
of km s _1 ), which could not detected at the SINFONI reso- 
lution (see Fig. 6). The velocity dispersions are rather high 
for a disk (compared to regular spirals), with V/a ~ 1.5-2 
in the observations and V/a ~ 2.5 in the model. The small 
difference is likely explained by the larger beam size in ob- 
servations. Models of galaxy mergers usually result in much 
lower V/a ratios, smaller than 0.5 in the case of multiple 
mergers (e.g., Bournaud et al. 2007). 

UDF 6462 is presently a quiescent galaxy. Its star forma- 
tion rate to stellar mass ratio SFR/M* ~ 1.5 Gyr -1 , which 
is consistent with the typical SFR — M* relation found for 
galaxies of this mass by Daddi et al. (2007) at z ~ 2. Stellar 
population studies in EE05 however suggest that the star 
formation activity was somewhat stronger in the last few 
10 8 yr. The formation of massive clumps in gas-rich disks 
can indeed trigger star formation: models in BEE07 typi- 
cally reach star formation rates of 30-40 Mq yr" 1 , up to 
ten times higher than the progenitor axisymmctric disks. 
Immeli (2004a) reports even higher star formation rates of 
200 Mq yr -1 but in the more extreme case of initially gas 
pure disks nearly as massive as the Milky Way. The life- 
time of these systems in BEE07 is about 0.5-1 Gyr, which 
is actually a lower limit because regulation by supcrnovae 
feedback was not included in these models. The circular 
velocity and dispersion estimated above corresponds to a 
dynamical (total) mass of ~ 4 — 5 x 10 10 Mq within the 
optical extent of the system, consistent with the stellar 
mass (Section 2.1) and implying a rather moderate (~20%) 
gas fraction at the observed instant 7 . However, an SFR of 
~ 50 Mq yr" 1 and an age of ~ 3 — 4 x 10 s yr imply an 



7 assuming that the mass within the optial radius is mostly 
visible with little contribution from the halo, as is generally the 
case for Local disk galaxies 
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Fig. 10. Position- Velocity diagram of UDF 6911 along its appar- 
ent major axis, in the [OIII] emission line. The thick pseudo-slit 
used to compute this diagram is shown in Fig. 2. 



initial gas fraction of about 50%, so that these mass esti- 
mates are compatible with the starting conditions of our 
model. We also note that this 50% initial gas fraction in 
the model is in agreement with observations by Erb et al. 
(2006b) and Daddi et al. (2008) and assumptions in Bouche 
et al. (2007). Furthermore, models with different initial gas 
fraction in BEE07 also show disturbed rotation comparable 
to the present case. 

Although the present models are initially asymmetric 
and highly complex structures develop quickly, real galax- 
ies could be even more complex, especially in the off-plane 
dimension, as a result of active inflow from cosmic streams 
and previous mergers. However, the model clearly indicates 
how internal dumpiness in a gas-rich galaxy can affect the 
kinematics. The comparison with observed properties then 
shows that the scenario in which the UDF 6462 clump clus- 
ter would have formed by internal instabilities in a gas-rich 
disk is a fully plausible one, in spite of the complex mor- 
phology and disturbed kinematics. These observed proper- 
ties can be explained, and the global rotating expected in 
this scenario is observed as a large-scale velocity gradient 
throughout the clump cluster. 



4.2. Kinematics of the bent-chain UDF 6911 

The Eastern companion of UDF 6462, UDF 6911 (Fig. 2), 
has the elongated morphology of a bent chain galaxy at 
z = 2.07. The same question on its origin can be posed as 
for UDF 6462. The kinematics of this system is only poorly 
resolved in our S1NFONI data. The P-V diagram of the 
[OIII] emission line along the apparent major axis of this 
object is shown in Fig. 10. Here again, we detect a veloc- 
ity gradient with a shape compatible with the flat rotation 
curve expected for a disk galaxy. Assuming that the circular 
velocity is traced by the envelope at 50% of the peak level 
(e.g., Donley et al. 2006) to account for beam-smearing ef- 
fects, and that this system is observed close to edge-on, the 
rotation velocity is estimated to be about 75 km s _1 . The 
resolution is too low to unveil local disturbances, but at 
least the global rotation pattern expected for a fragmented 
disk is found here, indicating that this interpretation is fully 
possible for this second object as well as for UDF 6462. 



4.3. Clump clusters and spiral disk and culge progenitors, 
and their role in the star formation history 

Within the plausible scenario where UDF 6462 and 
UDF 6911 (and possibly other clump clusters and chain 
galaxies) formed by internal fragmentation of clumpy disks, 
primordial galactic disks have to accrete enough mass at 
high redshift so that large fragments of their disks can 
collapse though gravitational instabilities. The presence of 
large gas reservoirs around disk galaxies at z ~ 1.5, un- 
veiled recently by Daddi et al. (2008), shows that large 
amounts of gas representing at least half of the total bary- 
onic mass can indeed be accumulated around primordial 
galactic disks. Massive disks with such high gas fractions 
are then likely to be gravitationally unstable and fragment 
into clump-cluster galaxies, or chain galaxies when seen 
edge-on, clump-clusters being face-on counterparts of the 
later as already suggested by Elmegreen et al. (2004). 

Fragmentation of gas-rich disks into clump clusters can 
then be an efficient mechanism to trigger star formation, 
with an efficiency at least comparable to that of major 
merger, if not higher, and longer star-formation timescales. 
Indeed, numerical models of merger-induced starbursts in- 
dicate that, on average, these have only modest intensities 
and durations (Di Matteo et al. 2007) even if some particu- 
lar cases can be very efficient (e.g., Cox et al. 2006). Recent 
observations further suggest that galaxy interactions and 
mergers have a modest influence on the cosmic star for- 
mation history (Jogee et al. 2008) and that internal pro- 
cesses in gas-rich disks play a more important role (Daddi 
et al. 2008). In a sample of Luminous Infrared Galaxies at 
0.4 < z < 1.2, Zheng et al. (2004) have identified a ma- 
jority of disk and irregular galaxies and a relatively small 
fraction of mergers. 

Later on, the clumps migrate towards the central region, 
being also progressively disrupted by the tidal field and 
their mutual interactions and redistributing the underlying 
disk mass into an exponential disk, as shown by numeri- 
cal simulations (Carollo et al. 2007, BEE07). Bulges with 
masses up to ^30% of the baryonic mass can be formed, or 
grown if already present, by the central merging of clump 
remnants; the bulges formed this way have the properties 
of the classical bulges of spiral galaxies (Elmegreen et al. 
2008b). Associated central black holes can be accreted at 
the same time from intermediate mass black holes formed 
by stellar collisions in the core of clumps, possibly explain- 
ing the black-hole-bulge scaling relations without mergers 
(see Johansson et al. 2008, for the case of mergers). Clump 
cluster and chain galaxies, when formed by the fragmen- 
tation of gas-rich primordial disks, are thus possible pro- 
genitors of present-day bright disk galaxies, as already sug- 
gested by van den Bergh et al. (1996). 

Clump clusters were found by BEE07 to form a thick 
disk of old stars, together which a thinner disk of gas and 
younger stars, which can further grow in mass by continued 
accretion. One can notice that the rotation velocity found 
in UDF 6462 is lower than what would be expected for the 
stellar mass from the Tully-Fisher relation (Conselice et al. 
2005). As the thick disks of spiral galaxies are observed to 
have a lack of rotation compared to thin disks (Yoachim 
& Dalcanton 2005, 2006), clumpy evolution at high red- 
shift is also a viable formation mechanism for thick disks, 
and indeed high redshift disks are observed to be thick 
(Reshetnikov et al. 2003, Elmegreen & Elmegreen 2006b). 
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5. Conclusion 

The UDF 6462 clump-cluster galaxy at z = 1.57 has a 
complex asymmetrical morphology, dominated by massive 
star-forming clumps. Spectroscopic observations of the Ha, 
[Nil] and [SII] emission lines of this system using SINFONI 
on UT4 at ESO/VLT have revealed a velocity gradient 
throughout the system, possibly tracing global rotation 
with a circular velocity of about 100 km s _1 . A central peak 
in the velocity dispersion further suggests that this veloc- 
ity gradient likely traces rotation. This is expected if this 
clump cluster formed by internal fragmentation of an unsta- 
ble gas-rich primordial disk. The velocity field shows large 
irregularities and dispersions which would be an anomaly 
for classical spiral disks, but can be explained by interac- 
tions between clumps, as shown by our numerical model. 
The unstable-disk model presented here reproduces, under 
a chosen projection, both the main morphological and kinc- 
matical properties of UDF 6462. 

An internal origin for clumpy disks is thus fully compat- 
ible with observed morphological and kinematical proper- 
ties in the case of UDF 6462. This does not definitely rule 
out a merger origin for this clump cluster. Nevertheless, 
the rotation-like velocity field would not be necessarily ex- 
pected this way, while it is predicted for internal instabil- 
ities in a gas-rich disks (BEE07). Furthermore, UDF 6462 
shows signs of central enrichment in a bulge-like central 
condensation, and a global metallicity in agreement with 
the estimated disk mass. This also is consistent with inter- 
nal dumpiness in a massive, gas-rich disk, and would be an 
unlikely coincidence for a merging system. A similar large- 
scale rotation-like velocity gradient is also found, although 
poorly resolved, in the nearby bent chain galaxy UDF 6911 
at redshift z = 2.07. 

Just like morphological asymmetries and dumpiness, 
large kinematical disturbances of several tens of km s _1 
can be a natural consequence of the internal evolution of 
gas-rich primordial disks if they become unstable and frag- 
ment into large star-forming clumps once they have ac- 
creted enough mass. Spectroscopic data should be obtained 
on more clumpy galaxies, and a systematic comparison with 
model predictions should be performed, before general con- 
clusions on the origin of chain and clump cluster galaxies 
can be obtained. At least, our present results show that 
even a severely disturbed case like UDF 6462 could actu- 
ally be formed by internal evolution and that hierarchical 
mergers are not necessarily needed to explain the forma- 
tion of such systems. Highly disturbed morphologies and/or 
disturbed kinematics with high velocity dispersions have of- 
ten been attributed to mergers in the recent literature. We 
have here shown that such properties do not provide ro- 
bust evidence for mergers, and may have an internal origin 
as well. Resolved spectroscopic data are needed to really 
constrain the origin of clumpy galaxies at high redshift on 
large samples, as well as large sets of numerical models 
to check that systematic kinematical classifications as in 
Flores et al. (2006) or Shapiro et al. (2008) can distinguish 
fragmented disks from merging systems. 

If clumpy galaxies are formed by large instabilities in 
gas-rich primordial disks, as is likely the case for UDF 6462 
and possibly also for UDF 6911, they can show large dis- 
turbances in their velocity field from one clump to the next 
one, but models presented here and in BEE07 predict that 
on the largest scales, velocity gradient tracing the initial ro- 
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Fig. A.l. Spectral energy distributions of the whole system (top 
line), central bulge- like object (middle line), and several clumps 
(bottom left). Symbols are the best model fits for the bulge 
and for one of the clumps. A good fit is obtained for the clumps 
without an old stellar population, while the bulge has a ~ 1 Gyr 
old population. The 'red wisp' is a small and somewhat redder 
clump at (a=3h 32m 43.68s - S=-27° 46' 58.8") close to the 
bulge itself. The typical uncertainties on the clump individual 
fluxes are 0.2 on \og(vf v ) in the optical and 0.3-0.4 in the NIR, 
so that the difference between the central bulge and the other 
clumps is robust. 



tation should be observed together with high velocity dis- 
persions. High turbulent speeds observed in high-redshift 
disks (Genzel et al. 2006; Forster Schreiber et al. 2006) are 
in agreement with this scenario, as they increase the Jeans 
mass and enable the formation of large, massive clumps in 
gas-rich disks. Moreover, UV and optical star forming disks 
at high-redshift have high angular momentum compared to 
local disks (Bouche et al. 2007), which phases of clumpy 
evolution can contribute to dissipate rapidly while fueling 
bulges and shaping more concentrated disks with exponen- 
tial profiles. The spectroscopic data presented in this paper 
thus support the hypothesis by Noguchi (1999), Immeli et 
al. (2004a,b), Elmegreen et al. (2005) and BEE07 that gas- 
rich primordial disks internally evolve through a clumpy 
phase into bright early-type disk galaxies with a massive 
exponential disk, a classical bulge, and possibly a central 
black hole. 



Appendix A: Stellar population study of UDF 6462 
bulge and star-forming clumps 

The observed SEDs for the whole UDF 6462 galaxy, the 
bulge, and several clumps are shown in figure A.l, using 
NICMOS resolution for the bulge and ACS resolution for 
the clumps. The bulge is clearly redder than the clumps. 
The reddest clump is an extended region jutting off from the 
bulge and may be part of the bulge. Theoretical fits to the 
SED of one of the clumps and to the bulge SED (crosses) 
are shown; details of these models are in EE05. Extinction 
curves from Calzetti et al. (2000) and Leitherer et al. (2002) 
were used, and intcrgalactic absorption by Lyman lines was 
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included following Madau (1995). The basic population evo- 
lution models are from Bruzual & Chariot (2003) assuming 
a metallicity of 0.008, which is 0.4 solar, and a Chabrier 
(2003) IMF. Star formation is assumed to begin at some 
time prior to the redshift of the object and to decay expo- 
nentially with a characteristic time scale. The age of the 
region and the decay time, along with the extinction, are 
parameters of the fit. The best fits here require extinctions 
of Ay = 0.6 mag for the clump and Ay = 1.2 mag for the 
bulge. The models suggest that the central, reddish, bulge- 
like object has an age of around 1 Gyr (or even older if 
the actual extinction is the same as for the clumps), and 
the clumps have an age of ~ 50 Myr with a comparable 
decay time in the star formation rate. There appears to 
be no significant old stellar component in the clumps. If 
these clumps were separate galaxies in a merging process, 
they would have older stellar components in them and look 
redder. No other part is as red as the central bulge or the 
reddish wisp near it (see also Fig. 3), which further sug- 
gests that no clump is a separate galaxy. The bulge is bluer 
than bulges in moderate redshift galaxies (e.g., Zheng et al. 
2004), which suggests some on-going star formation; how- 
ever the relatively poor resolution of NICMOS compared 
to the bulge clump observed with the ACS allows younger 
structure to contaminate the bulge SED, which would bias 
the bulge color slightly towards the blue. 

Appendix B: Metallicity distribution in galaxy 
mergers and fragmented disks 

The models from BEE07 used to interprete the kinematics 
of UDF 6462 do not directly provide predictions on the gas 
metallicity. To support the interpretation of the observed 
metallicity distribution in Fig. 7 and Section 4.1, we show 
some Tree-SPH models including chemical evolution from 
Di Matteo et al. (2007), for clumpy star- forming systems 
formed by galaxy mergers and internal disk fragmentation. 

B.l. Galaxy merger models 

The galaxy merger models shown here are from Di Matteo 
et al. (2007) and the Calmer 8 database (Chilingarian et 
al. in preparation). The merging disk galaxies start with a 
realistic radial metallicity gradient in their disk. These bi- 
nary mergers often have two main star forming regions at 
the center of each galaxy, but massive star forming regions 
also form in at larger radii by instabilities in spiral arms or 
in tidal tails (see e.g., Elmegreen et al. 1993), sometimes 
giving a clump-cluster morphology in the gas distribution. 
The four cases shown in Fig. B.l are early stage mergers 
that have been selected in the database for their clumpy 
morphology without any criterium related to the metallic- 
ity. They are representative of the majority of the ~ 900 
models available in the database. 

As shown in Fig. B.l, clump-cluster-like systems formed 
in early stages of galaxy mergers do not have a well or- 
ganized metallicity gradient. In UDF 6462, the metals are 
observed to be centrally concentrated around the red proto- 
bulge and star formation takes place in the outer clumps. At 
the opposite, in these galaxy mergers models the star form- 
ing knots in the outer regions are associated to metallicity 
peaks. Even the outermost large star-forming regions in 
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Fig. B.l. Gas density and metallicity maps in several galaxy 
merger models with a clumpy gas distribution. Outer star- 
forming regions can have been previously enriched in the inner 
disks of the progenitor spirals, like the tail to the bottom-left in 
case b, the clump to the left in case c, or the various external 
clumps in cases a and d. These cases are typical of the merger 
models from the Galmer database. Snapshots show between 
25 and 40 kpc in radius and are all in log scale, with different 
intervals for each one as they represent different systems. 



available at http://galmer.obspm.fr 



tidal tails around colliding galaxies, the TDG progenitors, 
are observationally known to have higher metallicities than 
classical disks outskirts, with quite often nearly-solar metal- 
licitities (e.g., Due et al. 2000). Indeed, star-forming regions 
at large radii in mergers can have been previously enriched 
in inner regions and expelled by tidal forces. Hence, the ra- 
dial metallicity gradient typical of disks is not preserved in 
the early phases of galaxy mergers (at least until the final 
relaxation) even when these have gas distribution similar 
to clumpy galaxies. 
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Fig. B.2. Clump-cluster model formed by internal disk frag- 
mentation. The gas density and metallicity maps are shown at 
two stages in the evolution. Green contours on the metal maps 
are metallicity contours, linearly spaced from to 0.05. Each 
box is 40 x 40 kpc. The star forming clumps have a moderate 
contribution to the metallicity distribution. The radial gradient 
from the progenitor disk is better preserved than in the merging 
of initially independent galaxies. 



B.2. Disk fragmentation models 

The model shown on Fig. B.2 is a gas-rich unstable disk, 
that gets fragmented in the same process as our dedicated 
model reproducing the kinematics of UDF 6462. Clumps 
form stars and contribute to the ISM enrichment, they in- 
teract and migrate inwards. Despite, a prominent metallic- 
ity peak remains clearly associated to the mass center, and 
no major metallicity peak is associated to the outer star 
forming clumps. The radial metallicity gradient typical of 
classical disks, and observed in UDF 6462, is hence best 
preserved if clump cluster galaxies form by the fragmen- 
tation of gas rich disks than in the merger of independent 
galaxies. 
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